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Abstract 

We present detailed experimental and numerical studies of random lasing in weakly scattering 
systems. The interference of scattered light, which is weak in the passive systems, is greatly 
enhanced in the presence of high gain, providing coherent and resonant feedback for lasing. The 
lasing modes are confined in the vicinity of the pumped volume due to absorption of emitted light 
outside it. In the ballistic regime where the size of gain volume is less than the scattering mean 
free path, lasing oscillation occurs along the direction in which the gain volume is most extended, 
producing directional laser output. The feedback for lasing originates mainly from backscattering 
of particles near the boundaries of pumped region. It results in nearly constant frequency spacing 
of lasing modes, which scales inversely with the maximum dimension of the gain volume. 

PACS numbers: 42.55.Zz,42.25.Dd 
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I. INTRODUCTION 



Random laser represents a non- conventional laser whose feedback is mediated by ran- 
dom fluctuation of dielectric constant in space. Since the pioneering work of Letokhov and 
coworkers^, lasing in disordered media has been a subject of intense theoretical and ex- 
perimental studies^. Random lasers have been realized in various material systems, from 
semiconductor nanoparticles, ceramic powder to polymers, organic materials and biologi- 
cal tissues. Their low fabrication cost, sample speciflc lasing frequency, small size, flexible 
shape, and substrate compatibility lead to many potential applications^. 

There are two kinds of feedback for random lasing: one is intensity or energy feedback, the 
other is field or amplitude feedback^i. The fleld feedback is phase sensitive (i.e. coherent), and 
therefore frequency dependent (i.e. resonant). The intensity feedback is phase insensitive 
(i.e. incoherent) and frequency independent (i.e. non- resonant). Based on the different 
feedback mechanisms, random lasers are classifled into two categories: (i) random laser 
with incoherent and non-resonant feedback, (ii) random laser with coherent and resonant 
feedback. 

In a strongly scattering system, multiple scattering facilitates light (of wavelength A) to 
return to the same coherence volume (~ A^) it has visited before, providing fleld feedback 
for lasing. The lasing frequencies are determined by the interference of scattered light 
returning via different paths. However, lasing with fleld feedback was realized also in the 
weak scattering regim o^i'^i^i^i'^'i^i^i'^" . In 1999, Frolov and coworkers reported laser-like emission 
from several weakly scattering samples including the vr-conjugated polymer fllms, organic 
dyes-doped gel fllms, opal crystals saturated with polymer and laser dye solutions under 
strip pumpin g^i^i^ . They observed two pump thresholds, the flrst corresponding to spectral 
narrowing of amplifled spontaneous emission (ASE), and the second to the transformation 
of featureless ASE spectrum into a flnely structured spectrum having numerous randomly 
positioned narrow lines. The flne spectral structures, which were reproducible under constant 
pumping condition and sample position, resulted from lasing with resonant feedback. The 
weak scattering, occurring repeatedly along the entire length of the excitation strip, supplied 
the distributed feedback for lasing. In 2004 Mujumdar and coworkers reported narrow 
emission spikes from suspensions of ZnO particles in Rhodmaine 6G-methanol solutions 
over a broad range of scattering strengths^. The spikes were distinct from shot to shot and 
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thus intrinsically stochastic. They were attributed to amplification of spontaneous emission 
along very long trajectories instead of optical cavities. Recently Poison and Vardeny used the 
power Fourier transform technique'' to reveal the underlying periodicity of the emission peaks 
in several disordered organic gain media in the weak scattering regimelS'. The ensemble- 
averaged power Fourier transform of random laser emission spectra contained a sharp, well- 
resolved Fourier component and its harmonics, which was characteristic of a well-defined 
laser resonator. However, it was not clear how such resonator was formed. 

The above puzzles motivated us to study random lasing in the weak scattering regime. 
The main questions we intend to address in this papers are: (i) is the interference effect still 
important to random lasing in the weak scattering regime? (ii) what are the lasing modes 
in the weakly scattering systems and how are they formed? To answer these questions, we 
performed both experimental studies and numerical simulations. 

The paper is organized as follows. A detailed experimental study of lasing in weakly 
scattering systems is presented in section II, followed by the numerical simulation in Section 
III. Section IV has a brief discussion and conclusion. 

II. EXPERIMENTS 

We performed experiments on several weakly-scattering systems which consisted of pas- 
sive scatterers embedded in active homogeneous media. The scatterers included Ti02 par- 
ticles of radius 200 nm, ZnO particles of radius 38 nm, and Si02 particles of radius 220 nm. 
The nanoparticles were suspended in a laser dye solution, e.g., Rhodamine 640 perchlorate, 
stilbene 420 or LDS 722 in diethylene glycol (DEG) or methanol. The experimental results 
obtained with different particles, dyes and solvents were qualitatively similar. As an exam- 
ple, we will demonstrate the lasing phenomena with colloidal suspensions of Ti02 particles 
in DEG with Rhodamine 640. 

A small amount of Ti02 (rutile) particles, with an average radius of 200 nm, were dissolved 
in the DEG solution of rhodamine 640 perchlorate dye. To prevent flocculation, the Ti02 
particles were coated with a thin layer of AI2O3. DEG was chosen as the solvent instead of 
the widely-used methanol because of the facts that (i) the windows of the quartz cuvette 
that contained the methanol solution were coated with a layer of Ti02 particles, whereas 
such coating was not observed for the DEG solution; (ii) the viscosity of DEG was about 
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30 times larger than methanol, thus the sedimentation of Ti02 particles in DEG was much 
slower. In our experiment, the particle density p ranged from 1.87 x 10^ cm"^ to 5.6 x 10^° 
cm~^. The scattering mean free path was estimated by Ig = i/pCs, where ag is the scattering 
cross section of a Ti02 spherical particle with radius 200 nm. The value of Ig varied from 
1.07 cm to 35 pm. The dye molarity M also changed from 3 to 10 mM. Right before 
the lasing experiment, the suspension was placed in an ultrasonic bath for 30 minutes to 
prevent sedimentation of the particles. During the experiment, the solution was contained 
in a quartz cuvette that was 1.0 cm long, 1.0 cm wide and 4.5 cm high. The dye molecules 
in the solution were optically pumped by the frequency-doubled output (Ap — 532 nm) of a 
mode-locked Nd:YAG laser (25 ps pulse width, 10 Hz repetition rate). The pump beam was 
focused by a lens into the solution through the front window of the cuvette. The radius of 
the pump spot at the entrance to the solution was about 20 /im. The experimental setup is 
shown schematically in the inset of Fig. 1. The emission from the solution was collected in 
the backward direction of the incident pump beam. A second lens focused the emission into 
a fiber bundle (FB) which was connected to the entrance slit of a spectrometer with cooled 
CCD array detector. The spectral resolution was 0.6 A. 

We started the experiment with a sample of M = 5 mM and p = 3.0 x 10^ cm~^. At 
low pumping level, the emission spectrum featured the broad spontaneous emission band of 
rhodamine 640 molecules. Above a threshold pump intensity, discrete narrow peaks emerged 
in the emission spectrum, and their intensities grew rapidly with increasing pumping. This 
behavior corresponded to the onset of lasing. The lasing peaks could be as narrow as 0.12 
nm. Their frequencies changed from pulse to pulse (shot). We repeated the experiment with 
samples of different particle density but the same dye concentration. Lasing was observed 
only within certain range of particle density. Figure 1 show the spectra of emission from 
five samples taken at the same incident pump pulse energy 0.4 pj. Each spectrum was 
integrated over 25 shots. Only a relatively broad amplified spontaneous emission (ASE) 
peak was observed for the neat dye solution, whereas a few discrete lasing peaks emerged 
on top of the ASE spectrum at small particle concentration p= 1.87 x 10^ cm~^. Increasing 
particle density to 1.87 x 10^ cm^"' led to an increase in the number of lasing peaks and the 
peak intensity. However, when p increased further to 1.3 x 10^° cm^'^, the lasing emission 
started to decrease. Eventually at p = 5.0 x 10^° cm~^ lasing peaks disappeared. A further 
increase of the incident pump pulse energy to 1.2 pj resulted in an ASE peak at longer 
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wavelength, shown in the right inset of Fig. 1. The red shift of the ASE peak might be 
caused by the surface effect on emission frequency of Rhodamine 640 molecules adsorbed on 
the Ti02 particles. One support for this explanation was that the emission frequency was 
blue shifted when we replaced the Ti02 particles by Si02 particles. 

In Fig. 2, the incident pump pulse energy at the lasing threshold Pt is plotted against 
the particle density p. At the lasing threshold, the slope of emission intensity versus pump 
pulse energy exhibited a sudden increase (inset of Fig. 2). At p = 3.8 x 10^ cm~^, lasing 
started at 0.21 //J. As p increased to 1.5 x 10^ cm~^, Pt decreased gradually to 0.12 //J. Then 
it remained nearly constant with a further increase of p. The threshold started to rise at 
p = 1.9 X 10^° cm^"^, then went up quickly with p. At p = 5.6 x 10^° cm^"^, no lasing peaks 
were observed up to the maximum pump pulse energy 2.0 pJ we used, although at 1.0 pJ 
an ASE peak appeared at longer wavelength. 

The particles played an essential role in the lasing process in our suspensions because 
lasing did not happen in the neat dye solution. One possibility was lasing within individual 
particles that served as laser resonators. It contradicted two experimental observations: (i) 
the lasing threshold depended on the particle density (Fig. 2); (ii) the laser output was 
highly directional (shown next). The left inset of Fig. 3 is a sketch of our directionality 
measurement setup. A fiber bundle was placed at the focal plane of the lens. It was scanned 
with fine steps parallel to the focal plane. At each step, the spectrum of emission into a 
particular direction was recorded. The output angle 9 was computed from the fiber bundle 
position. Its range was hmited by the diameter of the lens to about 14 degree. ^ = 
corresponded to the backward direction of the incident pump beam. In each spectrum, the 
emission intensity was integrated over the wavelength range of 604 — 612 nm in which the 
lasing peaks were located. Figure 3 is a plot of the integrated emission intensity versus the 
output angle 9 at p — 3.0 x 10^ cm~^. Although the spontaneous emission at low pumping 
was isotropic, the lasing emission was strongly confined to the backward direction of the 
incident pump beam. The divergence angle of the output laser beam was merely 4 
degree. To check the effect of refiection by the front window of the cuvette on lasing, we 
rotated the cuvette around the vertical axis and repeated the measurement. As shown in 
the right inset of Fig. 3, (f) represented the angle between the incident pump beam and 
the normal of the front window. Similar lasing phenomena were observed except a small 
increase of the lasing threshold. The lasing emission was always confined to the backward 
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direction of pump beam even when was much larger than the divergence angle of the 
focused pump beam, which was about 4 degrees. This result demonstrated that the front 
window of the cuvette was not indispensable to the lasing process. Figure 3 also shows the 
angular distribution of ASE from a sample of higher particle density (p = 5 x 10^° cm~^). 
The integrated intensity of ASE (at longer wavelength) was nearly constant over the angular 
range of detection. 

To understand the directionality of lasing emission from the dilute suspension, the 
pumped region was imaged through a side window of the cuvette. The measurement setup 
is sketched in the inset of Fig. 4(a). Emission from the excited region was collected through 
the side window by a 5x objective lens and imaged onto a CCD camera by integrating 
multiple pulses. Spectrum was taken simultaneously by partitioning the signal with a beam 
splitter (BS). Figure 4(a) compares the spectrum of emission through the side window to 
that through the front window of the cuvette from the same sample (p = 3 x 10^ cm~^, 
M — 5 mM) under identical pumping condition. The spectrum of emission from the front 
window exhibited large lasing peaks. However, only spontaneous emission was observed 
through the side window, and it shifted to longer wavelength as a result of reabsorption 
in the unpumped solution between the excited region and side window. To calibrate the 
reabsorption, we measured the spontaneous emission spectra at low pump intensity from 
both front and side windows. The magnitude of reabsorption was estimated from the in- 
tensity ratio of emission through the side window to that through the front. Based on this 
estimation, we concluded that the reabsorption was not strong enough to make the lasing 
peaks, which emerged in the front emission spectrum at high pump intensity, disappear in 
the side emission spectrum. This conclusion confirmed the result of lasing directionality 
measurement in Fig. 3. More importantly, the image of spontaneous emission intensity 
distribution taken through the side window exhibited the shape of the excited region in the 
sample. As shown in Fig. 4(b), the excited volume at low particle density had a cone shape. 
The length of the cone was much larger than its base diameter. Unfortunately, we could 
not get the exact length of the cone from the image, because near its end the spontaneous 
emission was too weak to be recorded by the CCD camera. At high particle density, the 
shape of excited volume changed to hemisphere as shown in Fig. 4(c) at a higher pump- 
ing power. This change was caused by increased scattering of pump light. In Fig. 4(b) 
p = 3.0 X 10^ cm"^, the scattering mean free path Ig at the pump wavelength Ap = 532 
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nm was estimated to be 800 /im. The (linear) absorption length la, obtained from the 
transmission measurement of neat dye solution, was about 50 nm. Strong 

pumping in the lasing experiment could saturate the absorption of dye molecules, leading 
to an increase of la- Since the shape of the excited volume shown in Fig. 4(b) was nearly 
identical to that in the neat dye solution, the scattering of pump light must be much weaker 
than absorption, i.e., la was still shorter than Is- In Fig. 4(c) p — 5.0 x 10^° cm~^, Ig was 
shortened to 53 /im. Scattering of pump light became much stronger. As a result of multiple 
scattering, the cone was replaced by a hemisphere. The image of excited volume provided 
some clue to high directionality of lasing emission at small p and non-directionality of ASE 
at large p in Fig. 3. At low particle density, stimulated emission in the cone-shaped gain 
volume was the strongest along the cone due to the longest path length. Since the cone 
was parallel to the incident pump beam, lasing was confined to the direction parallel to the 
pump beam. The divergence angle of laser output was determined by the aspect ratio 
of the excited cone, namely, ~ 2rp/ Lp, where Lp is the cone length and is the base 
radius. At large p, emitted photons experienced multiple scattering while being amplified 
in the hemisphere-shaped gain volume. Hence, the ASE was nearly isotropic. 

Therefore, the shape of gain volume determined the lasing directionality, i.e., lasing oc- 
curred along the direction in which the gain volume was most extended. However, it was 
still not clear how the laser cavities were formed in the dilute suspension. We examined the 
lasing spectra more carefully by taking single-shot emission spectra with the setup shown 
in the inset of Fig. 1. Surprisingly, in most single-shot spectra the spectral spacing of 
the lasing peaks was close to a constant. Figure 5(a) is an example of single-shot emission 
spectrum taken from the sample of p = 1.87 x 10^ cm~^ and M = 5 mM. The spectral 
correlation function C{dX) = (/(A)/(A + (iA))/(/(A)^) was computed for the spectrum in 
Fig. 5(a) and plotted in Fig. 5(b). The almost regularly spaced correlation peaks revealed 
the periodicity of lasing peaks. Despite that the lasing peaks completely changed from shot 
to shot, the peak spacing was nearly the same. In the spectrum taken over many shots, the 
periodicity was smeared out due to random (uncorrelated) peak positions in different shots. 
As an example, in Fig. 1 the solid thick curve represents the lasing spectrum integrated 
over 25 pulses, and the lasing peaks did not exhibit clear periodicity. We also noticed that 
the periodicity was less obvious at higher particle density. 

We would like to point out that the lasing peaks in the dilute suspension of particles are 
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fundamentally different from the stochastic ASE spikes that could be observed also in the 
neat dye solution. Figure 6 shows a single-shot spectrum of emission from the DEG solution 
of 5 mM Rhodamine 640 without any particles. The spectrum was taken under the same 
condition as that in Fig. 5. At high pumping level, stochastic spikes appeared on top of the 
ASE peak. The spikes in Fig. 6 were denser and narrower than the lasing peaks in Fig. 5. 
They changed constantly from shot to shot. When integrating the spectrum over subsequent 
shots, the spikes were quickly averaged out, leaving a smooth ASE spectrum shown as the 
dotted curve in Fig. 1. Note that the stochastic spikes also appeared in the spectrum of 
emission from the dye solutions with particles. However, they were taken over by the huge 
lasing peaks at high pumping level. The stochastic structure of the pulsed ASE spectrum 
was first reported thirty years agoM^. Since then, there have been detailed experimental 
and theoretical studies on this phenomeno n^^i^^i^^i^^ . The spectral fluctuation originated 
from random spontaneous emission, which was strongly amplified as it propagated through 
the cone-shaped pump volume. Non- stationary interference of the partially coherent ASE 
not only presented a grainy spatial pattern, but also caused drastic temporal fluctuation 
of intensity. The random intensity fluctuation within an ASE pulse in the time domain 
generated stochastic spikes in the spectral domain. According to the Fourier transformation, 
the width of spectral spikes was inversely proportional to the ASE pulse duration. In our 
case of picosecond pumping, the ASE pulse duration was of the order 25 ps. Thus the 
average spike width should be ~ 0.05 nm, which was close to the measured value of 0.07 
nm. 

To find out the location and size of the laser cavities in the dilute suspension, we placed 
a metallic rod in between the excited volume and the back window (inset of Fig. 7). On one 
hand, the rod prevented the emission from being reflected by the back window into the gain 
volume. On the other hand, a Fabry-Perot cavity was formed. Lasing in this cavity produced 
equally spaced peaks in the emission spectrum. The peak spacing AA was determined by 
the cavity length d (the distance between the rod and the front window), A A = \^ /2ned, 
where He is the effective index of suspension. As plotted in Fig. 7, AA decreased with 
increasing d. However, when d exceeded a critical value c/q — 450 /xm, AA jumped to a 
constant value and did not change with d any more (Fig. 7). The lasing spectrum and peak 
spacing aX d > d^ were identical to those from the same sample without the metallic rod. 
This result demonstrated that the laser cavity in the colloidal solution was located within 
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450 /im from the front window, i.e., in the vicinity of excited cone. Despite the scattering 
mean free path Ig (~ 800 fim) was much longer than the wavelength, the laser cavity was 
not extended over the entire sample, but confined to a region of dimension less than Ig. This 
could be explained by the reabsorption of laser emission in the unpumped part of sample. 
Our white-light absorption measurement and photoluminescence measurement of neat dye 
solution showed that Rhodamine 640 molecules in DEG had significant overlap between the 
absorption band and emission band. At the dye concentration M = 5 mM, the absorption 
length Ig at the emission wavelength Ag ~ 610 nm was about 300 fim. At low density of Ti02 
particles in the solution, the absorption length of emitted photons in the unpumped region 
was shorter than the scattering length. If the emitted light traveled beyond one absorption 
length from the pumped region, its chance of returning to the pumped region was extremely 
low. Therefore, the reabsorption of emission suppressed the feedback from the unpumped 
region of the system, and effectively reduced the system (or cavity) sizoii. 

The above experiment illustrated the laser cavity was confined in the vicinity of the 
excited region. The directionality of laser output suggested the laser cavity was oriented 
along the excited cone. Nearly constant spectral spacing of lasing modes resembled that of a 
Fabry-Perot cavity. Using the formula of a Fabry-Perot cavity, we derived the cavity length 
Lc from the average wavelength spacing AA of the lasing peaks, Lc = A^/2neAA ~ 200 — 300 
/im. The estimated cavity length was close to the length of excited cone observed from 
the side images. The cone length was determined by the penetration depth Lp of pump 
light into the suspension. When the scattering mean free path Ig was much larger than the 
absorption length at the pump wavelength, Lp was determined solely by absorption of pump 
light. This could explain the experimental observation that the averaging spacing of lasing 
peaks barely changed when the scattering mean free path Ig was varied by more than one 
order of magnitude. The cavity length Lc (or the penetration depth Lp) did not depend on 
Ig as long as Ig exceeded the absorption length. This result was confirmed by the side images 
of excited cones. 

In the case of linear absorption, Lp should be on the order of the (linear) absorption 
length la- However, the pumping in the lasing experiment was so intense that it saturated 
the absorption of dye molecules. The saturation photon flux density Ig = 1/afTf, where 
Cf is the fluorescence cross section and Tf is the lifetime of dye molecules in the excited 
state. For Rhodamine 640 molecules, o"/ is of the order 10~^^ cm^, and Tf 10~^ s. Thus, 
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Is ~ 10^^ cm~^s~^. The typical pump pulse energy at the lasing threshold was ~ 0.1 //J. 
Prom the pump pulse duration and pump spot size, we estimated that the incident pump 
photon flux density Ip ~ 10^^ cm-^s"-*^, which was two orders of magnitude higher than Ig. 
Hence, at the lasing threshold the absorption of dye molecules near the front window of the 
cuvette was strongly saturated, and the penetration depth Lp was much longer than the 
linear absorption length 

To conflrm Lc ~ Lp, we changed Lp by varying the dye concentration. The higher the 
dye concentration, the shorter the penetration depth. If L^. ~ Lp, the spacing of lasing peaks 
A A should increase. Figure 8(a) shows the single shot lasing spectra from three solutions of 
M = 3, 5, 10 mM. The particle density p was fixed at 3 x 10^ cm^^. It was evident that the 
spacing of lasing peaks increased at higher dye concentration. Images of the excited cones in 
the inset of Fig. 8(b) directly show that the excited cone was longer in the solution of lower 
dye concentration. Figure 8(b) plots the spectral correlation functions for the three spectra 
in Fig. 8(a). From them we extracted the average wavelength spacing AA = 0.34, 0.48, 0.96 
nm for M = 3, 5, 10 mM. This result conflrmed that the laser cavity length was determined 
by the pump penetration depth. Hence, the laser cavity was located in the excited cone, 
and the cavity length was approximately equal to the cone length. 

Although the lasing phenomenon in the weakly-scattering dye solution seemed to resem- 
ble that in a Fabry-Perot cavity with the "mirrors" at the tip and base of the excited cone, 
the question remained what made of the mirrors in the dilute colloidal suspension. The 
experiments presented earlier in this section already ruled out the front and back windows 
of the cuvette as the mirrors. The mirrors could not be formed by the nonlinear change of 
refractive index (both its real and imaginary parts) of the dye solution under intense pump- 
ing, otherwise lasing would have occurred also in the neat dye solution without particles. 
Another possibility is that the particles aggregated in the solution to form large clusters that 
served as mirrors. We monitored the solution during the experiment by imaging it onto a 
CCD camera through a side window of the cuvette. No clusters of size larger than l/im were 
observed in the pumped solution. We removed the clusters smaller than 1 /xm by filtering 
the solution. After the filtering, lasing phenomena remained the same. We also tried other 
suspensions such as ZnO in DEG and Si02 in methanol. Our previous studies confirmed 
the absence of particle aggregation in these suspensions, but we still observed similar lasing 
phenomenon. Therefore, the "mirrors" were not clusters of particles. Another candidate was 
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bubbles or shock waves that were generated by the pump pulse^^ii^^. We indeed observed 
bubbles in the solution when the pump beam was very strong and its focal spot was close 
to the front window of the cuvette. When the bubbles were big enough to be seen with our 
imaging apparatus, lasing peaks disappeared. Hence, big bubbles did not facilitate lasing. 
Small bubbles, which were invisible, might be generated when the pumping was not very 
high. Such bubbles were usually generated at the focal spot of the pump beam where the 
pump intensity was the highest. Thus the bubble formation should be sensitive to the dis- 
tance between the focal spot and the front window, which affected the pump intensity at the 
focal spot due to absorption in the solution. When we shifted the focal spot by moving the 
lens, the lasing behavior remained the same. This result eliminated the possibility of small 
bubbles contributing to lasing. All the experimental results led us to the conclusion that 
the coherent feedback for lasing came from the particles in the solution, more specifically, 
from the particles located near the tip and base of the excited cone. However, there were 
many particles inside the pumped volume, e.g., at p = 3 x 10^ cm~'^ the number of particles 
inside the excited cone was about 400. Why the feedback from the particles near the two 
ends of the cone dominated over that from the particles inside the cone? To understand this 
phenomena, we performed numerical calculations to be presented in the next section. 



III. NUMERICAL SIMULATION 



Several models were set up in the theoretical studies of stimulated emission in active 
random media, e.g. the diffusion equation with gain^LS^ and the Monte Carlo simulation^i^^. 
These models calculated light intensity instead of electromagnetic field, thus they ignored 
the interference effect. Although it is usually weak in the diffusive regime, the interference 
effect is not always negligible. One example is the coherent backscattering, namely, the 
interference between the counter-propagating light enhances the backscattered intensity by 
a factor of two. The experimental results in the last section illustrated that the interference 
effect was significant in weakly scattering samples with gain, leading to lasing with resonant 
feedback. In the dilute suspensions of particles in which lasing could be realized, the gain 
volume had a cone shape and the cone length was shorter than the scattering mean free 
path. If we considered only one photon, it most likely would not be scattered as it traveled 
from one end of the cone to the other. However, the intense pumping generated a huge 
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number of emitted photons. Despite the low probabihty of one photon being scattered, a 
significant number of emitted photons were scattered by the particles in the excited cone. 
Some of them were scattered backwards, providing feedback for lasing along the cone. Such 
weak feedback was greatly amplified as the backscattered light propagated along the cone. 
The interference of the backscattered light determined the lasing frequencies. 

Therefore, in the presence of large gain, the interference of scattered light is not negligible 
even in weakly scattering samples. To include the interference effect, we directly calculated 
the electromagnetic field in a random medium by solving the Maxwell's equations using the 
finite-difference time-domain (FDTD) method^^. The optical gain was modeled as negative 
conductance 



(Jo determined the gain magnitude, ti^o and I/T2 represented the center frequency and width 
of the gain spectrum, respectively. We neglected the gain saturation, and limited our calcu- 
lations to the regime just above the lasing threshold. A seed pulse with broad spectrum was 
launched at t = to initiate the amplification process. The lasing threshold was defined by 
the minimum gain coefficient (do) at which the electromagnetic field oscillation built up in 
time. 

Our numerical calculations aimed not at reproducing the experimental results, but at 
addressing the key issues and providing physical insight to the lasing mechanism. Hence, 
we simulated lasing in two-dimensional systems in order to shorten the computing time. To 
model the elongated gain volume in the experiment, optical gain was introduced to a strip 
of length Lp and width Wp. The refractive index was set at 1.0 both inside and outside 
the strip. Dielectric cylinders of radius 100 nm and refractive index 2.0 were introduced as 
scattering centers. The strip dimension was much smaller than the experimental value due 
to the limited computing power. Consequently, the number of cylinders inside the strip was 
reduced to keep the system in weak scattering regime. 

The experimental results in the last section suggested that the coherent feedback for 
lasing resulted mainly from the particles located near the tip and base of the excited cone. 
In the dilute suspension of particles, there was probably only one particle located at the tip 
of the cone. One question was whether the backscattering of a single particle could provide 
enough feedback for lasing. To answer this question, we started with only two particles in 
the gain strip, one at each end, in the numerical simulation. The total system size was 16 




I + i{uO - U0q)T2 I + i{uO + U0q)T2 



(1) 
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fim X 8 /im, and the gain strip 8 fim x 4 fim. At the system boundary, there was a perfectly 
matched absorbing layer. Two cylinders of radius 100 nm, placed at the two ends of the gain 
strip, had a separation of 8 /im. When the gain coefficient ctq was above a threshold value, we 
observed lasing oscillation. The spatial distribution of lasing intensity revealed that lasing 
occurred along the strip with the feedback from the two particles. The emission spectrum, 
obtained by Fourier transform of the electric field, consisted of multiple lasing modes equally 
spaced in frequency. We repeated the calculation after reducing the separation between the 
two particles inside the gain strip while keeping the strip length constant. The frequency 
spacing of lasing modes scaled inversely with the particle distance. These results confirmed 
lasing in the resonator composed of only two scatterers. 

In 1998, Wilhelmi proposed a laser composed of two Rayleigh scatterers with gain medium 
in between^. We generalized the Rayleigh scatterers to Mie scatterers and derived the lasing 
threshold condition: 



(Tf, is the backscattering cross section of one particle, which depends on the particle size, 
refractive index and light wavelength. The cavity length is the separation between the 
two scatterers. m is the dimensionality of the scattering system, is the threshold gain 
coefficient for lasing. ab/L"^'^ describes the probability of a photon being backscattered by 
one particle and propagating to the other particle. The larger the Lc, the less percentage 
of the backscattered photons can reach the other particle. It seems to suggest that the 
quality factor of the two-particle cavity decreases with increasing Lc- This perception is 
incorrect. At the lasing threshold, the cavity loss is equal to the gain, namely, the loss per 
unit length a = Qg. In the absence of intracavity absorption, a is related to the cold-cavity 
Q as a = 1/QX. From the threshold gain coefficient in Eq.Q, we derive the cold-cavity 
quality factor Q = Lc/{X\{m — l) InLc — lncrb|). As Lc —>■ oo, the numerator in the expression 
of Q diverges faster than the denominator, thus Q ^ oo. The rise of Q with Lc is attributed 
to the increase of one-path length of light inside the cavity. The reduction in the probability 
of photons backscattered by one particle then reaching the other particle in a long cavity is 
offset by the increase of photon travel time from one particle to the other. Hence, the lasing 
threshold decreases with increasing L^. The resonator with the lowest lasing threshold is 
composed of two particles with the largest possible separation inside the gain volume. 
Our numerical simulations also demonstrated directional lasing output from the two- 
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scatterer cavity. The near-field to far-field transformation of electric field gave the output 
laser intensity as a function of polar angle. Figure 9 shows the numerical data for three 
gain strips of length Lp = 4, 8, 16 fim. The strip width Wp was fixed at 4 yum. The two 
scatterers were always placed at the ends of the strip. From the envelop of the far-field 
intensity distribution, we obtained the angular width of the output laser beam. It decreased 
as the aspect ratio of the gain strip Lp/Wp increased. Similar results were obtained when we 
varied Wp and kept Lp constant. These results indicated that the output from a two-particle 
cavity laser cannot be simply regarded as scattering of a plane wave by a single particle even 
if Lc ^ A. It relied on both the geometry of the scatterers and the shape of the gain region. 
The directionality of the lasing output is a consequence of gain guiding. 

Although in the dilute suspension of particles there was probably only one particle at 
the tip of the excited cone, there were more particles at the cone base whose dimension 
exceeded the average distance between the particles. For example, when the particle density 
p = 5 X 10*^ cm~^ and pump spot radius ~ 20 fim, there were typically ten particles near the 
base of the excited cone. To simulate this situation, we placed ten scatterers randomly near 
one end of the gain strip and only one scatterer at the other end. The lasing peaks were 
almost equally spaced in frequency, with the spacing close to that with only two scatterers 
in the strip, one on either end. It suggested that the feedback from the ten scatterers near 
one end of the gain strip was equivalent to that from one located somewhere close to this 
end, as far as the lasing frequencies were concerned. 

One question we raised at the end of the previous section is why the feedback from the 
particles near the two ends of the excited cone dominated over that from the particles inside 
the cone. One possible explanation would be that light backscattered by the particles near 
one end of the cone experienced the most amplification as it traveled the longest path within 
the gain volume to the other end of the cone. If this were the reason, it implied the lasing 
modes would differ from the quasimodes of the passive system (without gain or loss). To 
check this conjecture, we randomly placed four cylinders inside the gain strip (8 /im x 4 
fim) in addition to the two at the ends, and compared the lasing modes to the quasimodes 
of the passive system. The quasimodes were calculated with the multipole method^^^^I. The 
field around each cylinder was expanded in a Fourier-Bessel series of regular and outgoing 
cylindrical harmonic functions. The Rayleigh identity related the regular part of the field at 
a particular cylinder to the waves sourced at all other scatterers. We found the quasimodes 
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by searching in the complex wavelength plane for the poles of the scattering operator. 

Our calculations revealed that despite the presence of additional scatterers, lasing still 
occurred in the direction parallel to the strip. Moreover, the lasing modes corresponded to 
the quasimodes with quality factor relatively high among all the quasimodes within the gain 
spectrum. Table 1 listed the frequencies v and quality factors Q of several lasing modes 
and the corresponding quasimodes. The gain spectrum was centered at 750 THz with a 
width of about 281 THz. The slight frequency shift of the lasing modes with respect to the 
quasimodes was due to the gain pulling effect. Figure 10 shows the intensity distributions 
of the lasing mode with v = 764.9 THz and the quasimode with u = 764.5 THz. It is 
evident that the lasing mode profile within the gain strip is nearly identical to that of the 
quasimode. We checked several lasing modes and obtained the same result. Therefore, the 
lasing modes were almost the same as the quasimodes in the presence of uniform gain. 

In Table 1, we also listed the frequencies and quality factors of the quasimodes in the 
two-particle case (without the four particles in the middle). Majority of the quasimodes in 
the system of six particles have frequencies similar to those of two particles. This comparison 
suggests in the weakly scattering system most quasimodes with relatively high Q are formed 
mainly by the feedback from the particles near the system boundary. The feedback from 
the particles in the interior of the system may slightly increase the quality factor or shift 
the mode frequency (see, e.g., the mode at u = 742.2 THz). However, the feedback from 
these particles may also be destructive and reduce the quality factor. As a result, the mode 
at z/ = 736.6 THz would not lase. Because most quasimodes with relatively high Q have 
frequencies similar to those with only the two particles furtherest apart, they, as well as the 
lasing modes, tend to be equally spaced in frequency. 

In the numerical simulation we did not place the particles outside the gain strip, because 
experimentally the feedback from those paticles is suppressed by reabsorption. Thus the 
effective system size is reduced, as shown in our previous calculatiooii. 

IV. DISCUSSION AND CONCLUSION 

Our experimental studies demonstrated lasing with field feedback in weakly scattering 
samples. The focused pump beam created a cone-shaped gain volume in the dye solution 
containing a small amount of nanoparticles. The cone length was determined by the ab- 
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sorption of dye molecules since optical scattering was much weaker than absorption. When 
the scattering mean free path exceeded the size of gain volume, lasing oscillation built up 
in the direction of strongest amplification, i.e., the direction in which the gain volume was 
most extended. This behavior was similar to that of amplified spontaneous emission in the 
weak scattering regime^^. The fundamental difference from ASE was, however, the exis- 
tence of feedback that originated from the backscattered light. As pointed out by Kumar 
and coworkers^S, the statistically rare sub-mean-free-path scattering could be made effective 
by strong amplification. In our experiment, the extreme weakness of feedback was compen- 
sated by high optical gain due to intense pumping. The interference of the backscattered 
light was greatly enhanced, leading to coherent and resonant feedback for lasing. 

It is important to note that the discrete lasing peaks were distinct from the stochastic ASE 
spikes. The latter originated from random spontaneous emission, which was significantly 
amplified in the presence of large gain. Such spectral fiuctuations also existed in ASE from 
the homogeneous media (without scattering). Hence, scattering was not indispensable to the 
existence of ASE spikes, although the fiuctuations could be enhanced by scattering which 
stretched the path length of photons inside the gain volume. In contrast, the presence of 
scatterers in the gain media was essential to the emergence of lasing peaks, indicating the 
lasing process relied on the feedback supplied by scattering. 

There have been many theoretical investigations of lasing modes in the weak scattering 
regime. They can be either extended modes^° or anomalously localized modea^. If only a 
small part of the random medium is pumped, the anomalously localized modes that locate 
inside the pumped region experience more gain than the extended modes that spread over 
the entire system. Therefore, the anomalously localized modes may have lower lasing thresh- 
old. However, the anomalously localized states are extremely rare in the weakly scattering 
samples. An alternative mechanism for spatial localization of lasing modes is absorption 
of emitted light outside the pumped region^i. The reabsorption of emission suppresses the 
feedback from the unpumped part of the random system and effectively reduces the system 
size. The lasing modes are therefore drastically different from the quasimode of the passive 
system (without gain or absorption). Even if all the quasimodes of the passive system are 
extended across the entire system, the lasing modes are still confined in the vicinity of the 
gain volume. 

The numerical simulations in section HI also illustrate that the lasing modes are nearly 
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identical to the quasimodes of the reduced system. The quasimodes are formed by dis- 
tributed feedback from all the particles inside the reduced system. The conventional dis- 
tributed feedback (DFB) lasers, made of periodic structures, operate either in the over- 
coupling regime or the under-coupling regime^. The random lasers, which can be regarded 
as randomly distributed feedback lasers, also have these two regimes of operation. In the 
under-coupling regime the quasimodes are formed mainly by the feedback from the scat- 
terers near the system boundary, while in the over-coupling regime the feedback from the 
scatterers inside the system becomes important. Thus, the quasimodes of a under-coupled 
system, especially the ones with relatively high quality, have almost regular frequency spac- 
ing. Note that the feedback from the scatterers inside the system is weak but not negligible, 
e.g., it may cause a slight shift of mode frequency or modification of the quality factor. 
In our experiment with the dilute suspensions of particles, the random lasing was in the 
under-coupling regime as a result of weak scattering and small size of the reduced system. 
Therefore, the dominant feedback from the particles near the cone ends resulted in nearly 
constant frequency spacing of the lasing peaks, which scaled inversely with the cone length. 
Due to weak feedback from the particles inside the cone, the lasing modes were not exactly 
equally spaced in frequency and some modes failed to lase as their quality factors were 
reduced. In the previous studies, e.g. Refs.i^*^*^^, the strong scattering or large system 
size make the random laser operate in the over-coupling regime. Thus the lasing peaks are 
randomly spaced. 

The results of our studies not only illustrated the physical mechanism of random las- 
ing in the weak scattering regime, but also demonstrated the possibility of controlling the 
frequencies and output directionality of random lasers by varying the pumping geometry, 
the scattering mean free path, and the absorption length at both excitation and emission 
wavelengths. Such control is important to the application of random lasers. 
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Table 1. Frequencies v (THz) and quality factors Q of lasing modes and quasimodes. 
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FIG. 1: Spectra of emission from DEG solutions of Rhodamine 640 (5 mM) and Ti02 particles. 
The particle densities are {l)p = (dotted curve), (2) 1.87 X 10^ cm-3 (thin solid curve), (3) 
1.87 X 10^ cm~"'(thick solid curve ), (4) 1.3 X 10^° cm-3 (dot-dot-dashed curve), (5) 5.0 x 10"'^'' cm *^ 
(dashed curve). All spectra were taken at the same pump pulse energy 0.4 /iJ. Each spectrum was 
integrated over 25 shots. Left inset is our experimental setup, BS: beam splitter, SP: spectrometer, 
FB: fiber bundle. Right inset is the emission spectrum shows the ASE peak, p = 5.0 x 10^'^ cm~'^. 
The pump pulse energy is 1.2 //J. 
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FIG. 2: The threshold pump pulse energy Pt as a function of the Ti02 particle densityp. The 
concentration of Rhodamine 640 in DEG is 5 mM. The inset is a plot of the emission intensity / 
versus the pump pulse energy Pe for the suspension with p = 3.0 x 10^ cm~^. 
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FIG. 3: The angular distribution of output emission intensity from DEG solutions of 5 mM Rho- 
damine 640 and 3.0 x 10^ cm~^ (solid square) or 5.0 x lO^'' cm~'^ (open circle) Ti02 particles. The 
solid line is a Gaussian fit of output laser beam. 9 = corresponds to the backward direction of 
pump beam. Left inset is a sketch of experimental setup, BS: beam splitter, SP: spectrometer, FB: 
fiber bundle. Right inset shows the angle (j) between the pump beam and the normal to the front 
window of the cuvette. 
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FIG. 4: (a): Spectra of emission through the side window (dashed curve) and front window (soUd 
curve) of the cuvette. The inset is a sketch of the experimental setup, BS: beam sphtter, SP: 
spectrometer, FB: fiber bundle, (b): Side image of excited region in DEG solution of 5 mM 
Rhodamine 640 and 3.0 x 10^ cm ^ Ti02 particles. The pump pulse energy is 0.2 /xJ. (c): Side 
image of excited region in DEG solution of 5 mM Rhodamine 640 and 5.0 x 10^° cm~^ Ti02 
particles. The pump pulse energy is 1.2 /^J. 
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FIG. 5: (a): A single shot emission spectrum from DEG solution of 5 mM Rhodamine 640 and 
1.87 X 10^ cm~^ Ti02 particles. The pump pulse energy is 0.2 /xJ pumping, (b): Spectral correlation 
function computed for the spectrum in (a). 
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FIG. 6: A single shot emission spectrum from the DEG solution of 5 mM Rhodamine 640 without 
any particles. The pump pulse energy is 0.3 fiJ pumping. 



26 




100 200 



300 400 500 

d (micron) 



600 700 



800 



FIG. 7: The wavelength spacing A A of lasing modes as a function of the distance d between the 
metalhc rod and the cuvette front window. The concentration of Rhodamine 640 in DEG is 5 mM, 
and the Ti02 particle density is 3.0 x 10^ cm^'^. The pump pulse energy is 0.2 /xJ. The inset is a 
sketch of experimental setup. BS: beam splitter, FB: fiber bundle. 
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FIG. 8: (a): Single shot emission spectra from colloidal solution of 3.0 x 10^ cm~^ Ti02 particles. 
The molarity of Rhodamine 640 in DEG is (from top to bottom) 3 mM, 5 mM, and 10 mM. (b) 
Spectral correlation functions of the single shot emission spectra in (a). Insets are side images of 
pumped region for in each solution. 



28 




-J I I I I I I I t i^-l I I I I I L. 

-80 -60 -40 -20 20 40 60 80 

Angle (degree) 



FIG. 9: Far-field intensity of laser emission as a function of polar angle. The zero degree corresponds 
to the direction parallel to the gain strip. The width of gain strip Wp was fixed at 4 /im. The strip 
length Lp = 4 /im (dotted curve), 8 /im (dashed curve), and 16 /zm (solid curve). The inset is a 
sketch of the geometry of the 2D system in our numerical simulation. Two cylinders of radius 100 
nm are located on both ends of the gain strip. 




FIG. 10: (a) Spatial intensity distribution of a lasing mode at = 764.9 THz in the 2D system of 
six dielectric cylinders inside a gain strip 8 /xm x 4 /xm. The circles represent the cylinders. The 
rectangle marks the boundary of the gain strip, (b) Spatial intensity distribution of quasimode at 
V = 764.5 THz in the same system as (a) but without gain. 
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